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Summary 

Bat species around the world have recently been recognized as major reservoirs of 
several zoonotic viruses, such as severe acute respiratory syndrome coronavirus 
(SARS-CoV), Middle East respiratory syndrome coronavirus (MERS-CoV), 
Nipah virus and Elendra virus. In this study, consensus primer-based reverse tran¬ 
scriptase polymerase chain reactions (RT-PCRs) and high-throughput sequencing 
were performed to investigate viruses in bat faecal samples collected at 11 natural 
bat habitat sites from July to December 2015 in Korea. Diverse coronaviruses were 
first detected in Korean bat faeces, including alphacoronaviruses, SARS-CoV-like 
and MERS-CoV-like betacoronaviruses. In addition, we identified a novel bat 
rotavirus belonging to group H rotavirus which has only been described in 
human and pigs until now. Therefore, our results suggest the need for continuing 
surveillance and additional virological studies in domestic bat. 
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Introduction 

Bats are considered reservoirs of several emerging infec¬ 
tious disease. Recently emerging human viruses, such as 
severe acute respiratory syndrome coronavirus (SARS- 
CoV), Middle East respiratory syndrome coronavirus 
(MERS-CoV), Nipah virus, Hendra virus and Ebola virus, 
are thought be bat-borne viruses (Han et al., 2015). In 
addition, new species of influenza A viruses, lyssaviruses, 
paramyxoviruses and coronaviruses continue to be discov¬ 
ered in bats around the world (Tong et al., 2013; Banyard 
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et al., 2014; Yang et al., 2014; Mortlock et al., 2015). Reas- 
sortant group A rotavirus was also detected in the faeces of 
a straw-coloured fruit bat (Esona et al., 2010). Therefore, 
viruses that originate in bats may be a source of additional 
spill over from wildlife into domestic animals and humans 
(Plowright et al., 2015). 

Since the outbreak of MERS-CoV in Korea in 2015, the 
prediction and control of newly emerging viruses has 
gained urgency. However, in Korea, there have been a few 
studies on viruses in domestic bat species. Serological evi¬ 
dence of Japanese encephalitis virus (JEV) was reported in 
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the early 1990s (Lee and Lee, 1992). In 1995, there was a 
report on the isolation and genetic characterization on a 
hantavirus from bats (Jung and Kim, 1995). More informa¬ 
tion of viruses in bats is needed to clarify the role of these 
animals in infectious diseases epidemiology. 

Therefore, in this study, we investigated viruses in bat 
species in Korea, using 49 faecal samples collected from July 
to December 2015 in 11 sites in natural bat habitats. 

Materials and Methods 

Samples 

From July to December in 2015, 49 faecal samples were col¬ 
lected at 11 sites in natural bat habitats, such as caves, an 
abandoned mine and under a bridge (Table 1). Fresh bat 
faeces were placed into transport medium in a 10% suspen¬ 
sion and were transported to the laboratory for further lab¬ 
oratory analysis. 

The major bat species at the collection sites were deter¬ 
mined based on morphology and on previous data from 
bats’ roosting sites (Han et al., 2012). In this study, the 
major bat species in the collecting sites were Rhinolophus 
ferrumequinum, Myotis macrodactylus, Myotis aurascens 
Kuzyakin, Myotis petax and Miniopterus schreibersii. 


Viral detection 

Consensus primer-based reverse transcriptase polymerase 
chain reactions (RT-PCRs) were performed to detect 
influenza A viruses, coronaviruses, lyssaviruses and param¬ 
yxoviruses (Fouchier et al., 2000; Poon et al., 2005; 
Vazquez-Moron et al., 2006; Tong et al., 2008). The 
detailed information on the RT-PCRs is in Table SI. A total 
of 49 bat faecal samples were tested. The amplicons from 
positive RT-PCRs were sequenced using target-specific for¬ 
ward and reverse primers synthesized by Cosmogenetech 
Co. Ltd (Daejeon, Korea). The nucleotide data obtained 
from Sanger sequencing of the PCR fragments were further 
analysed with related sequences in GenBank using BioEdit 
(Hall, 1999) and MEGA version 6 (Tamura et al., 2013). 

High-throughput sequencing 

In total, 16 bat faecal samples that were collected in four 
sites were used for high-throughput sequencing (Table 1). 
The supernatants from suspensions of one, four, five and 
six faecal samples collected in Danyang, Wonju, Taebaek 
and Uljin, respectively, were pooled by sampling site. Each 
of the four pooled samples was filtered through a 0.2-/(m 


Table 1 . Information on bat faecal samples 


Cases 

Collected date 

Region 

Site name 

Major bat species 

No. of samples 

1 

24 July 2015 

Danyang 

Under HSA Bridge 

Unknown 

6 C 

2 

4 August 2015 

Wonju 

HO cave 

Rhinolophus ferrumequinum, 
Myotis macrodactylus 

4 C 

3 

14 August 2015 

Taebaek 

YY cave 3 

Myotis aurascens Kuzyakin 

Myotis petax 

Rhinolophus ferrumequinum 

5 C 

4 

16 August 2015 

Uljin 

SL cave b 

Unknown 

1 c 

5 

9 October 2015 

Jeongsun 

PG cave 

Rhinolophus ferrumequinum 

4 

6 

14 October 2015 

Inje 

BM abandoned 

mine 

Rhinolophus ferrumequinum 

6 

7 

18 October 2015 

Taebaek 

YY cave 3 

Myotis aurascens Kuzyakin 

Myotis petax 

Rhinolophus ferrumequinum 

3 

8 

25 October 2015 

Bonghwa 

DRN cave 

Rhinolophus ferrumequinum 
Myotis macrodactylus 

3 

9 

18 November 2015 

Taebaek 

YY cave 3 

Myotis aurascens Kuzyakin 

Myotis petax 

Rhinolophus ferrumequinum 

1 

10 

19 November 2015 

Uljin 

SL cave b 

Unknown 

1 

11 

1 December 2015 

Hapcheon 

BT cave 

Unknown 

3 

12 

1 December 2015 

Munkyung 

BGS cave 

Miniopterus schreibersii 

4 

13 

1 December 2015 

Yeongwol 

YJ cave 

Rhinolophus ferrumequinum 

4 

14 

1 December 2015 

Mi won 

CS cave 

Rhinolophus ferrumequinum 

4 


Collected at the same site (YY cave) on a different date. 
b collected at the same site (SL cave) on a different date. 
c pooled and sequenced by high-throughput sequencing. 
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filter and ultracentrifuged at 252 000 g for 1 h. Each pellet 
was resuspended in 500 /d of lx digestion buffer (Turbo 
DNA Free Kit; Ambion, Darmstadt, Germany), and a 
total of 10 /d of Turbo DNase was added. After incubat¬ 
ing for 30 min at 37°C, suspensions were transferred to 
a 1.5-ml reaction tube and centrifuged at 1500 g for 
3 min. The supernatants were used for RNA extraction. 
RNA was extracted using Trizol LS (Invitrogen Corp., 
Carlsbad, CA, USA), following the manufacturer’s man¬ 
ual. The extracted RNA was submitted to Macrogen 
(Seoul, Korea) for high-throughput sequencing in a 
HiSeq 2000 sequencing system based on the transcrip- 
tome de novo sequencing platform. 

The obtained viral contigs were analysed and annotated 
by the MG-RAST server (Meyer et al., 2008). The cut-off 
for the annotation was 10~ 5 , 60% and 15 for maximum 
e-value, minimum percentage identity and minimum align¬ 
ment length, respectively. The contigs annotated as mam¬ 
malian viruses were validated through BLASTn (http:// 
www.ncbi.nlm.nih.gov). The GenBank accession numbers 
for the nucleotide sequences in this study are KU528584- 
KU528594. 

Results 

PCR-based detection of influenza A viruses, 
coronaviruses, lyssaviruses and paramyxoviruses in bat 
faeces 

There were no samples positive for influenza A viruses, lys¬ 
saviruses or paramyxoviruses. 

Among 14 cases of the submitted samples, three cases 
were positive for coronavirus by RT-PCR, showing the 
expected 440 nt target. These included the following: (Case 
no. 1) one of four samples from HO cave in Wonju where 
Rhinolophus ferrumequinum and Myotis macrodactylus were 
the major species present, (Case no. 6) one of six samples 
BM abandoned mine in Inje where Rhinolophus ferrume¬ 
quinum was the major species found and (Case no. 12) two 
of four samples BGS (Case no. 12) cave in Munkyung 
where Miniopterus schreibersii was found to be the major 
species (Table 1). 

In the HO cave, one of four faecal samples (B15-8) was 
positive by coronavirus-specific RT-PCR. When the PCR- 
positive fragment was sequenced, it was found to be a bat 
coronavirus (Bat CoV) belonging to the genus alphacoron- 
avirus and showing 89% nucleotide sequence identity with 
Bat CoV A620 (GenBank accession no. DQ648828) 
detected in China in 2005. 

In the BM abandoned mine, one of four faecal samples 
(B15-21) was found to be a Bat CoV belonging to the genus 
betacoronavirus and showing 98% nucleotide identity with 
SARS-like Bat CoV Lingbao-331 (GenBank accession no. 
KF294456) detected in China in 2012. 
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In the BGS cave, two of four faecal samples (B15-40 and 
B15-41) were found to be Bat CoVs belonging to the genus 
alphacoronavirus and showing 99% nucleotide identity 
with Bat CoV A895 (GenBank accession no. DQ648841) 
and 98% nucleotide identity with Bat CoV A910 (GenBank 
accession no. DQ648849), both of which were detected in 
China in 2005. 

High-throughput sequencing from the pooled bat faecal 
samples 

As shown in Table 1, cases 1, 2, 3 and 4 were further anal¬ 
ysed by high-throughput sequencing. By Hiseq 2000 
sequencing, 4341, 162 559, 7690 and 16 657 contigs were 
obtained from cases 1, 2, 3 and 4, respectively. Based on the 
annotation by the MG-RAST system, several viral 
sequences that are related to mammalian viruses, insect 
viruses, plant viruses, fungal viruses and bacteriophages 
were identified. 

The mammalian viral sequences annotated by the MG- 
RAST system were further validated by BLASTn (http:// 
www.ncbi.nlm.nih.gov). The BLASTn result is presented 
in Table 2. Fifteen viral sequences from case no. 1 were 
similar to porcine rotavirus SKA-1 (69-73% nucleotide 
identity), Bat CoV SC2013 (70-91% nucleotide identity), 
Bat CoV HKU5-5 (82% nucleotide identity), Alphacoron¬ 
avirus Eptesicus fuscus-related strains (95-96% nucleotide 
identity) and Rhinolophus pearsoni bunyavirus Shaanxi 
2011 (72% nucleotide identity). Three viral sequences 
from case no. 2 were similar to Banna virus 02VN078b 
(90% nucleotide identity), Bat CoV CDPHE15 (80% 
nucleotide identity) and Bat CoV Neixiang-31 (82% 
nucleotide identity). There was no mammalian viral 
sequences in cases 3 and 4. 

Genetic characteristics of bat coronaviruses in Korea 

Partial RNA-dependent RNA polymerase (RDRP) 
sequences obtained from B15-8 (sample in case no. 2), 
B15-21 (sample in case no. 6), B15-40 and B15-41 (sample 
in case no. 12) were phylogenetically compared with other 
reference sequences (Fig. la). Three of these sequences 
clustered with the alphacoronaviruses (B15-8, 40 and 41), 
and the other one clustered with the betacoronaviruses 
(B15-21). Bat CoV B15-8 grouped with bat CoVs found in 
China, showing 85.7-88.8% nucleotide identities, and was 
found to be closely related with Bat CoV-512 and porcine 
epidemic diarrhoea virus (PEDV)-CV777, showing 78.6% 
and 74.6% nucleotide identities, respectively. Bat CoV B15- 
21 grouped with SARS-CoV-like Bat CoVs, such as Bat 
CoV Rfl and Bat CoV HKU3 with 90.3-97.7% nucleotide 
identities, and was found to be closely related to human 
SARS-CoVs (SARS-CoV Tor2, MA15, ExoNl strains), 
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Table 2. Mammalian viral sequences matched 
BLASTn search parameters _ Matched strain with contigs from high-throughput sequencing 

Size % Nucleotide 


Contigs 3 

(bases) 

f-value 

identity 

Accession No. 

Strain name 

C6091 

404 

2e-37 

73 

AB576634 

Porcine rotavirus SKA-1 

C6719 

438 

1 e-21 

69 

AB576631 

Porcine rotavirus SKA-1 

C7554 

211 

9e-19 

73 

AB576625 

Porcine rotavirus SKA-1 

C9980 

215 

5e-22 

73 

AB576629 

Porcine rotavirus SKA-1 

Cl 0008 

285 

2e-35 

73 

AB576629 

Porcine rotavirus SKA-1 

Cl 0298 

355 

2e-48 

73 

AB576631 

Porcine rotavirus SKA-1 

C3339 

230 

1e-55 

83 

KJ473821 

BtVs-BetaCoV/SC2013 

C6171 

203 

1 e-16 

70 

KJ473821 

BtVs-BetaCoV/SC2013 

Cl 0151 

219 

1e-36 

82 

EF065512 

Bat coronavirus PIKU5-5 

Cl 0579 

207 

3e-67 

91 

KJ473821 

BtVs-BetaCoV/SC2013 

C6073 

243 

6e-98 

95 

HQ585081 

Alphacoronavirus Eptesicus 
fuscus/Appalachian Ridge/ 

PI-Cl 148/IT/USA/2009 

C8820 

417 

0.0 

96 

JX537911 

Alphacoronavirus Eptesicus fuscus 
strain NeCoV isolate Ffrag_ 
SAGA_bigbrown_2010 

C8271 

270 

2e-61 

95 

HQ585082 

Alphacoronavirus Eptesicus 
fuscus/Appalachian Ridge/ 

PI -C837/IT/USA/2009 

C7729 

209 

2e-21 

76 

KF294382 

Anlong Ms bat coronavirus 
isolate Neixiang-32 

C5770 

313 

2e-34 

72 

KC154063 

Rhinolophus pearsoni bunyavirus 
isolate Bat Rp-BunyaV/ 
Shaanxi2011 

C700585 

207 

8e-70 

90 

EU265697 

Banna virus isolate 02VN078b 

C817114 

201 

6e-40 

80 

KF430219 

Bat coronavirus CDPHE1 5/USA/ 

2006 

C695165 

226 

6e-53 

82 

KF294381 

Neixiang Md bat coronavirus 
isolate Neixiang-31 

3 Contigs obtained from case no. 1 and case no 2. No mammalian viral sequences were detected in 
case no. 3 and case no 4. 


showing 89.1% nucleotide identity. Bat CoV B15-40 and 
B15-41 grouped with Bat CoV 1A and IB, and nucleotide 
identities within the group were 92.3-97.1%. 

Because an RDRP sequence from betacoronavirus- 
related Bat CoV B15-1-6 (detected from the pooled sample 
in case no. 1) was not obtained, a partial spike gene 
sequence (cl0151) was compared with the spike genes of 
reference betacoronaviruses (Fig. lb). In addition, the full 
spike gene sequence obtained from Bat CoV B15-21 was 
also analysed. Bat CoV B15-1-6 grouped with MERS-CoV- 
like Bat CoVs (Bat CoV strains HKU5, A434, SC2013 and 
HKU4), showing 80.9-82.2% nucleotide identities, and was 
found to be closely related to MERS-CoV EMC, showing 
76.6% of nucleotide identity. 

The amino acids sequence encoded by the full spike gene 
of B15-21 was aligned and compared with human SARS- 
CoV and SARS-CoV-like Bat CoV sequences (Figure SI). 
The S2 domain is more conserved than the SI domain in 
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the spike protein. The receptor binding domain (RBD) of 
the SARS-CoV-like Bat CoVs, except for RaSITC014, had 
two major deletion sites, TGNYN (446-450) and 
PFSPDGKPCTPPA (472-484), compared to human SARS 
coronaviruses sequences. The receptor binding motif 
(RBM) and human ACE2 (hACE2)-interacting residue 
were highly variable between human SARS-CoVs and 
SARS-CoV-like Bat CoVs but RaSEIC014 strain. 

Phylogenetic characterization of Group H bat rotavirus in 
Korea 

Porcine rotavirus-related sequences from pooled bat faeces 
in case no. 1 were analysed with reference sequences of 
rotaviruses in human and animals (Fig. 2). Based on the 
partial sequences of the VP1, VP3 and VP4 regions, the 
rotavirus-related sequences were found to belong to rota¬ 
virus group H. The group H bat rotavirus in this study 
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(b) 


Alpha CoVs 


Gamma CoVs 



Bat CoV-HKU5(KC522101) 

Bat CoV-A434(DQ648790) 

Bat CoV-SC2 013(KJ473821) 
B15-1-6 | 

MERS-CoV(NC_019843) 

Bat CoV-HKU4-1(NC009019) 

Bat CoV-RsSHC014 (KC881005) 
SARS CoV-Smo 1-11 (AY485277) 
SARS CoV-Ufbam (AY278741) 
Bat CoV-Rp3 (DQ071615) 
B15-21 | 

Bat CoV-Rf1(OQ412042) 

Bat CoV-A701(OQ648793) 

HCoV-NL63(AY567487) 

HCoV-229E(AF304460) 

PEDV-CV777(AF353511) 

BatCoVA515(DQS48791) 

BatCoV-HKU 10(JQ989270) 


MERS CoV-like 
Clade 


SARS CoV-like 
Clade 


Fig. 1. Phylogenetic trees constructed from bat coronavirus sequences detected in this study and from other coronaviruses. (a) bootstrap consensus 
tree inferred from 1000 replicates with the UPGMA method (evolutionary distance parameter: maximum composite likelihood method) based on 351 
bases of the RNA-dependent RNA polymerase gene, (b) bootstrap consensus tree inferred from 1000 replicates using the neighbour-joining method 
(evolutionary distance parameter: maximum composite likelihood method) based on 163 bases of the spike gene. The box indicates the bat CoV 
sequences obtained in this study. 


showed 69.4-72.7%, 71.1-72.9% and 67.3-72.2% nucleo¬ 
tide identities with the other group H rotaviruses, based on 
VP1, VP3 and VP4 partial sequences, respectively. 

Discussion 

According to a recent report, there are 23 species of bats in 
Korea (Han et al., 2012). Most Korean bats are thought to 
be insectivores, as no fruit bats have been found in Korea 
to date, except for some fruit bats as pets or indoor exhibi¬ 
tions (Han et al., 2010). Through the target-specific RT- 
PCR, partial RDRP sequences of three alphacoronaviruses 
(B15-8, 40 and 41) and one betacoronavirus (B15-21) were 
detected in the faeces of Korean bats. Bat CoV B15-8 was 
detected in faeces collected at the HO cave where Rhinolo¬ 
phus ferrumequinum and Myotis macrodactylus were the 
major species. This virus was found to be closely related 
with Bat CoVs recently found in China and PEDV, a por¬ 
cine diarrhoeic virus around world. On the other hand. Bat 
CoVs B15-40 and 41 were detected at the BGS cave, where 
Miniopterus schreibersii was the major species. These 
sequences were closely related with Bat CoV 1A and IB, 
which were reported to have a clear host species restriction 
to Miniopterus magnate and Miniopterus pusillus in Hong 
Kong (Chu et al., 2008). As Bat CoVs B15-40 and 41 were 
also found in faeces in a cave where Miniopterus schreibersii 

© 2016 Blackwell Verlag GmbH • Transboundary and Emerging Diseases. 


was dominant. Bat CoV strains in the same lineage with 
Bat CoV 1A and IB may be adapted to the Miniopterus spp. 

It was interesting that Bat CoV B15-21 found in this 
study was closely related to SARS-CoV-like Bat CoVs, such 
as Bat CoV Rfl, based on the partial RDRP and S genes, 
respectively. The newly identified Bat CoV B15-21 clustered 
with the SARS-CoV-like Bat CoV group in the phylogenetic 
analysis. The major bat species in the site where Bat CoV 
B15-21 was collected was Rhinolophus ferrumequinum, 
which is consistent with a previous report that SARS-CoV- 
like bat CoVs were mostly found in Rhinolophus spp. (Yuan 
et al., 2010). The amino acids of its full spike gene showed 
that the B15-21 strain is closely related to SARS-CoV-like 
Bat CoVs Rp3, Rfl, HKU3-2 and 273 rather than human 
SARS-CoVs and Bat CoV RaSHC014, which was reported 
to have high potential for human emergence (Menachery 
et al., 2015) and showed two specific deletion sites in the 
RBD of the spike protein. One additional betacoronavirus 
(B15-1-6), which was not detected by the consensus pri¬ 
mer-based RDRP-specific RT-PCR (Poon et al., 2005), was 
detected by high-throughput sequencing as shown in 
Table 2. Three contigs were similar to Bat CoV SC2013, 
which is a MERS-related betacoronavirus that was found in 
bats in China (Yang et al., 2014). One contig was also simi¬ 
lar to Bat CoV HKU5-5, which clustered with MERS-CoV 
(Woo et al., 2007; Annan et al., 2013). Phylogenetic 
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B219(Human)-EF453355 
03V0567(Chicken)-NC_021590 
CAL-1 (Human)-EU490414 
Bang373(Human)-NC_021541 
03V0568(Chicken)-NC_021625 
05V0049(Chicken)-NC_014511 
KU(Human)-AB022765 
OSU(Porcine)-GU199514 
Bristol(Human)-NC_007547 
Cowden(Porcine)-M74216 
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RVB 
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J19(Human)-DQ113900 
B219(Human)-EF453357 
SKA-1 (Porcine)-AB 576631 
| Bat RVH (C10298)| 
03V0567(Chicken)-NC_021589 
Cal-1 (Human )-EU490417 
Bang-373(Human)-NC_021551 
KU(Human)-AB022767 
OSU(Porcine)-AY277921 
Bristol(Human)-NC_007574 
Cowden(Porcine)-M74219 
05V0049(Chicken)-NC_014514 
03V0568(Chicken)-NC_021630 
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RVG 
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RVA 
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RVD/RVF 


— J19(Human) 

— B219(Human) 

-[ Bat RVH(c7554)| 

— BR59(Porcine) 
-SKA-1 (Porcine) 

— 03V0567(Chicken) 

— IDIR(Murine) 

—ADRV(Human) 

— Cal-1 (Human) 

— Bang-373(Human) 

— KU(Human) 

— OSU(Porcine) 
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— Cowden(Porcine) 
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RVB 

RVA 

RVC 

RVD/RVF 


Fig. 2. Phylogenetic trees constructed from Group H rotavirus-related 
sequences detected in bat faeces in this study and sequences of other, 
representative rotavirus strains. The bootstrap consensus tree inferred 
from 1000 replicates using the neighbour-joining method (evolutionary 
distances were computed using the Kimura 2-parameter method), (a) 
275 bases of the VP1 region, (b) 347 bases of the VP3 region, (c) 202 
bases of the VP4 region. The box indicates the group H bat rotavirus 
sequences obtained in this study. 


analysis based on partial spike gene sequences of Bat CoV 
B15-1-6 showed that it belonged to a MERS-CoV-like clade 
with previously known MERS-CoV and MERS-CoV-like 
Bat CoVs. 

Although information of the associated bat species was 
not clearly available, this study first report that SARS-CoV 
and MERS-CoV-like bat CoVs exist in Korea. This may 
imply that other strains or types of those coronaviruses 
may circulate in bat species in Korea. Therefore, we cannot 
ignore appearance of novel or variant coronavirus through 
local emergence in Korea, like SARS and MERS in China 
and Middle East, respectively. 

In addition to coronaviruses, this study is the first to 
determine the existence of group H rotaviruses in bat faeces 
collected in Korea. Group H rotavirus is a recently pro¬ 
posed group of rotaviruses that include strains ADRV-N 


and B219, which infect human adults (Alam et al., 2007; 
Matthijnssens et al., 2012). So far, group H rotaviruses 
have only been reported in human and pigs (Molinari 
et al., 2015), but this study provides evidence that bat spe¬ 
cies may be a host of group H RVs. To confirm that, there 
should be follow-up studies including virus isolation and 
characterization, genomic analysis, continuous surveillance 
and VP6-based classification (Matthijnssens et al., 2012) to 
find its prevalence, epidemiology and zoonotic potential. 

Localized emergences of zoonotic diseases are usually 
due to expansion of the wildlife-human interface (Morse 
et al., 2012). The human-bat interface in Korea has not 
been fully described. People usually recognize that bats are 
flying and roosting in the deep forest or caves. Elowever, 
the sites where the bat faecal samples were collected showed 
some evidence of human-bat contact (Figure S2). Some 
caves or abandoned mine where bats are roosting are 
located near human habitats. People living there frequently 
visit into the caves or the abandoned mine to take a rest or 
to ferment foods. Some caves were also used for some reli¬ 
gious ceremony by Korean traditional shamans. 

In this study, SARS-CoV-like and MERS-CoV-like bat 
CoVs and group H rotavirus were detected for this first 
time in Korea, which may be of interest because of their 
zoonosis potential. Several kinds of evidence for human- 
bat contact were also found. Therefore, the potential for 
novel or variant human viruses from bats in Korea should 
not be ignored. Continuous surveillance and additional 
virological research should be a priority to prevent future 
emergences of zoonotic viral diseases. 
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